Interaction between ionic lattices and superconducting condensates 
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The interaction of the ionic lattice with the superconducting condensate is treated in terms of 
the electrostatic force in superconductors. It is shown that this force is similar but not identical to 
the force suggested by the volume difference of the normal and superconducting states. The BCS 
theory shows larger deviations than the two-fluid model. 
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I. INTRODUCTION 

The theory of deformable superconductors deals either 
with effects of the lattice deformations on the super- 
conducting condensate or with deformations of the crys- 
tal lattice driven by the inhomogeneous superconducting 
condensate. For example, a lattice deformation around a 
dislocation pins a vortex Ja^. In contrary, forces generated 
by supercurrents contribute to the magnetostriction^^ 
and a condensate depletion at the vortex core deforms 
the lattice so strongly that a significant renormalization 
of the vortex mass has been predicted jSiLSal In some 
cases one cannot say which of the two effects are domi- 
nant. This happens, for instance, if the ionic lattice de- 
formation influences the structure or orientation of the 
Abrikosov vortex latticed 

The free energy describing deformable superconduc- 
tors has to include at least three parts. The first part 
is the elastic energy of deformations. Its structure and 
parametrization have well been established for a long 
time^ The second part is the magnetic energy and the 
energy of superconducting condensation. This part can 
be covered on different levels. Here we will refer to the 
Ginzburg-Landau (GL) theory^ employed in the major- 
ity of the above mentioned studies. The third part is a 
cross term which describes the mutual effect of deforma- 
tions and the condensate. In this paper we focus on this 
interaction term. 

In the phenomcnological approach put forward by 
Kramer and Bauer- 1 the interaction of the lattice and 
the condensate is described by a local product of the 
lattice density with the density of superconducting elec- 
trons. If the lattice is modeled by an isotropic deformable 
medium and the interaction is assumed to be local, the 
interaction energy of Kramer and Bauer is the only one 
compatible with the system symmetry. Indeed, the con- 
densate density is a scalar which can interact only with 
another scalar. The only isotropic scalar quantity linear 
in the deformation is the trace of the strain tensor, which 
is proportional to the ionic density. Of course, one can 
construct more elaborate scalars within non-linear terms 



but these are higher-order corrections. 

Two modifications are at hand. First, one can take into 
account that the real lattice is never isotropic. Even in 
simplest lattices of elementary metals the shear rigidity 
depends on the orientation of the deformation with re- 
spect to the crystal axes. In the anisotropic crystal, there 
are two scalar quantities linear in the shear deformation. 
A corresponding anisotropic generalization of the inter- 
action between the condensate and the ionic lattice has 
been discussed by Kogan et al~& 

Second, one can go beyond the local approximation. 
The reason for such a step is the following. The in- 
teraction in the local approximation is justified only 
for neutral systems^ The superconducting condensate, 
however, drives the system out of neutrality induc- 
ing the electrostatic potential known as the Bernoulli 
potential, 1 ^ The non-local interaction mediated by the 
Bernoulli potential in the bulk of the superconductor has 
been discussed in our previous paper fi 5 - 

Another contribution to the charge transfer induced 
by the condensate is the surface dipoleJ^ While all the 
above mentioned interactions result in a force density 
acting in the bulk of the crystal, the surface dipole yields 
the force which acts as an external pressure imposed on 
the surface. As far as we know, the surface dipole has 
never been discussed within the theory of deformable su- 
perconductors. In this paper we want to fill this gap. 

The paper is organized as follows. In section II we 
show that the surface dipole determines changes of the 
crystal volume during its transition from the normal to 
the superconducting state. To this end we first introduce 
the basic concept in Sec. II. A and derive the coefficient 
of the local interaction from the pressure dependence of 
the condensation energy at zero temperature in Sec. II. B. 
The result is compared with the force due to the surface 
dipole in Sec. II. C. In section III we introduce the in- 
teraction mediated by the Bernoulli potential. We first 
derive a formula for the coefficient of the local interac- 
tion. In Sees. III. A and III.B we evaluate the interaction 
coefficient for moderately strong and weak coupling su- 
perconductors. In section IV we discuss differences and 
conclude. 
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II. LOCAL APPROACH 

In their pioneering study Kramer and Bauer— proposed 
to deduce the interaction strength from the pressure de- 
pendence of the critical magnetic field B c . Since exper- 
imental data for this parametrization are conveniently 
found in literature, this approach has been employed by 
other authors too. 

In this section we provide a derivation of the local in- 
teraction of Kramer and Bauer within the GL picture of 
the superconductor. The presented approach is based on 
papers by Simanek^ and Hake?^ 



A. Phenomenological force 

The volume of a metal changes at the phase transition 
from V n in the normal state to V s in the superconducting 
state i 17 ! 18 This change is described by a relative change 
a of the specific volume defined as 



V n — V s = aV 3 . 



(1) 



If the specific volume becomes inhomogeneous, the 
crystal has regions requiring different distances of neigh- 
boring atoms. This leads to internal stresses which can 
be expressed via an effective force density^ 

F ph = KVa, (2) 

where K is the modulus of hydrostatic compression or 
simply the bulk modulus. It is defined as the inverse of 
the relative volume change with respect to the pressure 

----- (3) 

k ~ v d P ■ [ } 

The temperature dependence of a is similar to the tem- 
perature dependence of the superconducting fraction 

(4) 

where subscript zero denotes the zero temperature value. 
Simanek&i and Coffey^ use the BCS gap A, 

while other authors prefer the GL function ip. We will 
restrict our attention to the vicinity of the critical tem- 
perature, where both forms are equivalent since the BCS 
gap and the GL function are linearly proportional to each 
other ^ 

Relation ([4J is the central approximation in the phe- 
nomenological theory of deformable superconductors. 
Assuming that in the normal state the system is ho- 
mogeneous, ao and ipo are constants. With the GL 
function normalized to the density of pair-able electrons 



2\H 



we obtain 



F ph = -a KV\4>\ 



(6) 



which we will use in our discussion. 



B. Difference of normal and superconducting 
volume 



Now we link the relative change of the specific volume 
a to the pressure dependence of the condensation energy 
e con . We follow the derivation of Hake»il 

The volume of the sample is the pressure derivative at 
fixed temperature of the Gibbs free energy 



dp 



(7) 



At zero magnetic field and zero temperature, the free 
energy of the normal state G n is higher than the super- 
conducting free energy G s by the condensation energy 



G n — G s — V s s c 



(8) 



Due to the complete expulsion of the magnetic field 
from type-I superconductors, the condensation energy is 
conveniently observed via the critical magnetic field at 
zero temperature 



2/j.q' 



(9) 



In his study, Hake expresses all thermodynamical rela- 
tions exclusively in terms of the critical magnetic field 
Here we prefer to use the condensation energy. 

^From the pressure derivative of the relation (§| it fol- 
lows 

8f 



v n -v s = v s - 



dp 



Comparing the thermodynamical relation (|I0|) with the 
definition (|TJ) we obtain the coefficient ao in terms of the 
condensation energy 



a 



de c 



1 dV s 



dp —v s dp- 

In terms of the bulk modulus ([3|) we have 



a 



de c 



dp 



K 



(11) 



(12) 



The force ([6]) depends on the product otoK, it is thus 
advantageous to introduce the inverse bulk modulus also 
into the first term of Eq. (|12p . For simplicity we consider 
hydrostatic pressure and conventional superconductors 
with isotropic structure. In this case we can express the 
pressure dependence of the density of the condensation 
energy e con via its dependence on the electron density, 



ds c 



dp 



de con dn 
dn dp 



(13) 



Since the number of electrons iV does not change, we can 
express the bulk compressibility via the change of the 
density n — N/V as 



1 

K 



1 dn 

n dp 



(14) 
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Using relation (fl~3"| and the bulk modulus (fT4|) in equation 
(H2) we get 



a K = n— e c 

on 



(15) 



The density derivative is taken under the condition of 
charge neutrality, i.e., the lattice density changes with 
the electron density. 

The phenomenological force ([6]) according to relation 
([12]) thus reads 



(16) 



We will compare this form with the electrostatic force 
resulting from a surface dipole. 



C. Compression via the surface dipole 

At the surface of a metal the electrostatic potential 
rises by few Volts from its vacuum value to the value 
deep in the metali 2 ^ This increase is spread partly out- 
side the region occupied by ions, typically on the scale 
of the tunneling length of electrons in the potential bar- 
rier given by the work function. A part of the barrier 
is located inside the metal on the scale of the Thomas- 
Fermi screening length. Both scales are of the order of 
Angstr0ms making the potential step very sharp. This 
sharp step is called the surface dipole. 

The surface dipole naturally depends on the tempera- 
ture. Moreover, when the metal undergoes a transition to 
the superconducting state, the temperature dependence 
of the surface dipole changes. Briefly, the superconduct- 
ing condensate affects the surface dipoleJ^ Let us denote 
this additional potential near the surface as tpx- 

An amplitude <Pt(0) — (pxipo) of the additional po- 
tential step follows from the Budd-Vannimenus theorem 
cole 



which changes the volume of the crystal by 

a V = -^-Pr- 
op 



(20) 



Clearly, the surface dipole contributes to the relative 
change of the specific volume a. From ([20)) we find 



a = -^piat [<pr(0) - </?t(oo) 



(21) 



At zero temperature / = — e con , therefore from (|17|) it 
follows 



Plat [<A>(0) ^ <A)(°o)] = n 



di 



(22) 



The dipole-induced volume change at zero temperature 
thus reads 



a 



1 / de c 



K \ dr 



- Sc 



(23) 



Comparing (|23|) with (|15|) one can see that the volume 
change is fully induced by the surface dipole 



Q = OtQ. 



(24) 



The fact that the volume change is driven by the sur- 
face dipole shows that one should be cautious using the 
relative change of the specific volume ag as a coefficient 
of the interaction between the ionic lattice and the con- 
densate. 

Studies of the electrostatic potential in superconduc- 
tors have shown that the bulk and surface potentials are 
of different nature and are covered by distinct theories. 
We note that these theories are experimentally verified. 
The surface potential including the surface dipole has 
been observed by Morris and Brown via the Kelvin ca- 
pacitive pickup i 14 ' 16 The internal charge transfer caused 
by the bulk electrostatic potential in the vortex core has 
been observed by Kumagai, Nozaki and Matsuda via the 
nuclear magnetic resonance ! 21 ' 22 



Plat [<Pt(0) - <pr(oo)] = / 



dn 



(17) 



where / = / s — /„ is the free-energy density by which the 
superconducting state differs from the normal state and 
where pi at is the charge density of the ionic lattice. We 
assume a superconductor which fills the half-space x > 0. 

The additional potential exerts an electrostatic force 
density on the ionic lattice 



Ft = -piatV^T- 



(18) 



The integral of this force density across the surface region 
corresponds to an effective pressure on the lattice 

oo 

Pt = J dxF% = Put [tp T {0) - <Pt(oo)] , (19) 



III. ELECTROSTATIC FORCE ON IONS 

According to the Hellmann-Feynman theorem, elec- 
trons act on ions exclusively via the electrostatic force. 30 
In this spirit, we expect the force density to be of elec- 
trostatic nature, 



(25) 



where ip is the electrostatic potential created by the su- 
perconducting electrons which is conveniently derived fol- 
lowing the approach of Rickayzen* 2 ^ Since the system is 
in equilibrium, the Gibbs electrochemical potential jj, for 
electrons is constant all over the sample. It is locally 
defined from the density of free energy / as 



fx = eip + 



df 
dn' 



(26) 
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Following the customary choice in the theory of super- 
conductivity we set the electrochemical potential to zero, 
/i = 0, therefore 



<P = 



1&1 

e dn 



(27) 



The theory of the electrostatic potential has been de- 
rived under the assumption that the ion lattice is stiff 
and its deformation is not included. A combination of 
both effects has not been studied so far, therefore it is 
not yet clear how the density derivative in (|27|) is mod- 
ified by lattice deformations. For simplicity we assume 
that the density derivatives in (f2"T)) and in (fT7|) are the 
same. This is the case if the pairing interaction has a 
purely electronic nature so that the lattice density has 
no effect on the condensation energy. 

Formula (|2"T|) is quite general. It has been employed 
by Rickayzen to evaluate the Bernoulli potential in super- 
conductors using the London theory supplemented by the 
phenomenological temperature dependence of the super- 
conducting density. In the same paper- 2 -^ Rickayzen has 
used formula (j27|) with the BCS free energy and recov- 
ered the result of Adkins and Waldram^ 

We use the free-energy density in the GL approxima- 
tion 



o(T-T c )M a + i&M 4 



2m* 



-ihV 



e *A)V| 2 + -^|B a - V x A| 2 , 

ZfJ,Q 



(28) 



where A is the vector potential and B a is the applied 
magnetic field. The reader not familiar with the GL the- 
ory is referred to the textbook 2 ^ of Tinkham. 

The magnetic free energy (the last term of (1281) ) does 
not depend on the electron density. For simplicity we 
also assume that the Cooper pair mass m* is indepen- 
dent of this density. Since the GL wave function ip and 
the vector potential A are independent variational fields, 
the density derivative of the free energy yields the elec- 
trostatic potential 



adT c 
e dn 



IVf 



T -T c da 
e dn 



1 db 

2e dn 



H 4 - (29) 



The nonlocal corrections discussed in Ref . UM are hid- 
den in the second and third terms of (f2"5)) . They can be 
made explicit using relations for material parameters a 
and b, e.g. (|33|) . and the GL equation, which couples 
nonlocal and nonlinear contributions. 

Here we restrict our attention to a close vicinity of the 
critical temperature, where all nonlocal and nonlinear 
contributions can be neglected. Indeed, for T — > T c the 
GL wave function vanishes |V>| 2 oc T c — T. In lowest order 
in T — T c we can neglect the second and the quartic term 
so that relation ([29]) simplifies to 



adT c 
e dn 



h/f 



Now we are ready to evaluate the force acting on the 
ionic lattice. The electrostatic force density (|25|) with the 
electrostatic potential (|30|) reads 



Piat^ / dT c 2 



(31) 



In first order in T c — T, gradients of material parameters 
a and dT c /dn do not contribute, i.e., 



(32) 



We have used p\ a t = —en demanded by the local charge 
neutrality. 

With nonlocal and nonlinear corrections neglected, the 
electrostatic force (|32p like the phenomenological force 
(|16p are proportional to the gradient of the supercon- 
ducting density \tp\ ■ Our next aim is to compare the 
electrostatic coefficient an^j- with its phenomenological 
precursor 2( d£ g° n — ^f 1 )- We will show that the relative 
values of these coefficients depend on the strength of the 
pairing interaction. 



A. Superconductors with moderately strong 
coupling 

To be able to compare the electrostatic force density 
with the phenomenological force density ^ , we need 
the GL parameter a as function of the electron density n. 
For metals like niobium or lead, it is possible to use the 
asymptotic form of the two-fluid free energy of Gorter 
and Casimir— giving^ 



ace 



(33) 



Here 7 is the linear coefficient of the specific heat. 

The critical temperature T c and the critical magnetic 
field Bo at zero temperature are linked via the conden- 
sation energy. The two-fluid model yields 2 ^ 



4 7 c 2u 



(34) 



Within this approximation, the electrostatic force density 
(1321) reads 



FGC 
el 



de c 



di 



4 c dn 



v|Vf 



(35) 



One can see that this is similar but not identical to the 
phenomenological force density Q16p . In particular, the 
same dominant term oc de COTL /dn results from both ap- 
proaches. 

We note that Simanek and many of other authors use 
the approximation | d£ Q° l n \ ^> l^f^j, i- e -, they consider 
only the derivative in their formulas. Within this accu- 



(30) 

racy both formulas are equivalent. 
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B. Superconductors with weak coupling 

Metals like aluminum have weak electron-phonon cou- 
pling and one can use BCS relations. This approximation 
results in a slightly different electrostatic force. 

/From the BCS theory Gor'kov has obtained parame- 
ters of the GL theoryJ^ The linear GL coefficient reads 



obcs 



67r 2 fc B T c 
7C(3)£f ' 



(36) 



where Ef is the Fermi energy. The Riemann Zeta func- 
tion has the value £(3) = 1.202. 

The BCS and the Gorter and Casimir approxima- 
tions of a can be related within the free electron model. 
The electron density determines the Fermi vector k-p — 
(37r 2 n) 1//3 in terms of which Ef = h 2 k F /2m. We also use 
7 = (2/3)7r 2 fc|iVo, where N = (l/4n 2 ){2m/h 2 )k F is the 
single-spin density of states. Combining these relations 
one finds that both values differ by a numerical factor 



obcs 



12 
703) 



«gc = 1-43 ace- 



(37) 



The BCS relation connecting the critical temperature 
with the condensation energy defined via the magnetic 
field yields another numerical factor 



4' c 2 Mo 

In the two-fluid model we find T C GC = 
correction 



(38) 



The 



T, 



BCS 



T C GC 



= 1.028 T, 



GO 



(39) 



is by an order of magnitude less important than the factor 
1.43 from relation (|37|) . however. 

/From (|38[) we obtain T c in terms of B , which we use 
in the force density (J35J. With the BCS relation ([go]), 
the electrostatic force density (|3"2"]l reads 



? BCS 



24 



1 de c 



7C(3) V 0.947 
de. 



On 



= 3.012 



dn 



-0.!)47^^Vl'1 2 ^"» 



One can see that the dominant contribution is increased 
by slightly more than 50% as compared to the phe- 
nomenological force density (|16[) and the Gorter-Casimir 
approximation (|35p . Compared to the Gorter-Casimir 
approximation the critical temperature T c is replaced 



additionally by the BCS value T C BCS . It will thus be inter- 
esting to test the validity of the phenomenological force 
on materials of rather different coupling strength. 
IV. DISCUSSION 



We have shown that the change of the volume dur- 
ing the superconducting transition can be expressed as 
a compression caused by the surface dipole. In analo- 
gy we have used the internal electrostatic potential for 
the lattice deformation in the bulk. The force resulting 
from the electrostatic potential in the superconductor is 
similar but not identical to the phenomenological force 
suggested from the volume change. 

For moderately coupled materials well described by 
the Gorter-Casimir two-fluid model, the phenomenolog- 
ical and the electrostatic forces have identical dominant 
terms oc ds con /dn. They differ in the correction terms 
only. For weakly coupled superconductors covered by the 
BCS theory, the dominant term is increased by nearly 
50%. We have not discussed the strongly coupled super- 
conductors which do not obey any of these limits. One 
can expect that the dominant term is reduced. 

It is a question whether the above derived small dif- 
ferences in the internal forces can be accessed by some 
of recent experiments. Among the physical phenomena 
mentioned in the introduction, the magnetostriction of- 
fers the most sensitive experimental technique. Indeed, 
one can resolve even deformations driven by such small 
changes in magnetization as those caused by the de Haas- 
van Alphen effect with relative changes of the suscepti- 
bility of the order of one in ten millions^ A supercon- 
ductor in the Meissner state is an ideal diamagnet with 
large but fixed magnetization. Small deviations from ide- 
ality appear due to the penetration of the magnetic field 
into the surface. At the same time, the magnetostric- 
tion combines the internal forces with the surface dipole. 
To separate these two contributions, it will be necessary 
to analyze how the deformation depends on the sample 
geometry. 

The two approaches compared in this paper represent 
two extreme models. In the electrostatic approach all 
forces are attributed to a mean electric field. In the phe- 
nomenological model the system is treated as locally neu- 
tral which implies that all forces are attributed to bonds 
between ions. We expect that a realistic description re- 
quires to combine both approaches. 
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